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Abstract 

It is proposed that the electron-pairing medium of the iron-based 
superconductors may be the orbital fluctuation of the transition metal ions. But the 
characteristic frequency of the orbital fluctuation has not been given. For the first time, 
the author has calculated the real-time evolution of the electron clouds of transition 
metal ions in BaFe2As2 under excitations by the time-dependent density functional 
theory (TDDFT). There are different modes of fluctuations. The characteristic 
frequencies are 150 meV, 160 meV, 250 meV, and 200 meV, respectively, for the 
modes the author observed. The results are unexpected, because the general view is 
that the change of the electron density is very quick, and the frequency is much higher 
than the lattice vibration. The frequencies the author obtained are close to that of the 
lattice vibration in conventional superconductors at normal and high pressures, 
indicating the orbital (or electron cloud) fluctuation can by the electron pairing 
medium. Based on the calculation results, the author proposed a new electron pairing 
mechanism. 
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1 Introduction 
The electron-lattice interaction cannot explain the electron-pairing mechanism in 


the unconventional copper-oxide |": 7! and iron-based P 4! superconductors. As a new 
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electron pairing mechanism, orbital fluctuation has been extensively studied ISH78], 
but the characteristic frequency of the orbital fluctuation has been not obtained. If the 
frequency is too low or too high, it is unlikely to be an electron-pairing medium. Table 
1 gives the maximum phonon frequencies and Te of three typical conventional 
superconductors. It is reasonable to guess that the characteristic frequency should be 
about 100-250 meV. 


Table 1 Maximum phonon frequencies and T. of three typical conventional superconductors. 


Pb [9] MgB> [9] H;S [10] 


a / meV 9 90 250 


T./K 7 39 164 


Based on the above consideration, the author has made an investigation, 
exploring the characteristic frequency of the orbital fluctuation. The real-time 
evolution of charge densities under excitations of BaFe2As2 has been calculated by 


TDDFT "!: !?] method. This paper will report the methods and results. 


2 Methods 


13, 14 


Calculations were conducted within the Octopus package |’? '4!, The generalized 


gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) was used to 


[15, 16] 


describe the exchange-correlation energy. The GGA+U method was considered 


to deal with the strong correlation of the 3d electrons. HSCV pseudopotentials H” 
were adopted. Approximated Enforced Time-Reversal Symmetry (AETRS) algorithm 
was used to approximate the evolution operator and the time step is 0.002 ii/eV. 
3p°d°4s? of Fe, 3d!°4s’p* of As and 5p°6s? of Ba are taken as valence electrons in the 
calculation. The real space spacing is 0.35 Bohr. The vibration of charge densities is 
excited with two methods. The first is a time-dependent electric field along 111 
direction: 


t—2t)- 


{=F cos(7 o coset If |t—t,|>t). then f(t) =0. 


0 


The second is to change some atom slightly. In the TDDFT calculation, the ions 


are fixed and the electron densities vibrate freely. 


2 


3 Results and Discussions 

Electric fields with two frequencies were used to excite the charge fluctuation. 
When œ = 1.8 eV (the applied electric field), Figure 1(a) shows the time dependence 
of the electric filed. The excitation vanishes after 5000 steps and the charge densities 
vibrate freely thereafter. After 8000 steps’ evolution, the charge densities are recorded. 
Figure 1(b) shows the crystal structure and the real-time evolution of the change of 
the charge density. 100 represents the charge density difference between the 8100th 
step and the 8000th step and the same below. The change of the electron clouds of Fe 
ions (marked by 3 and 4) become obvious gradually with the evolution steps. To give 
a better illustration, Figure 1(c) shows the maximum value of the charge density 
change as a function of the evolution steps. The charge density change reaches its 
maximum value after about 850 steps, and the corresponding time is 1.7 A/eV. It is the 
time from zero to the maximum. The time of one period should be 41.7 fi/eV and the 
frequency is about 150 meV. This result is unexpected, because the general view is 
that the change of the electron density is very quick, and the frequency is much higher 
than the lattice vibration. But the results obtained by the author are close to the lattice 
vibration. The frequency given in this way is not very accurate, but the frequency is 
consistent with the authors’ guess, indicating that it can be the electron-pairing 


medium in unconventional superconductors. 
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Figure l(a) Time dependence of the electric field applied, œ = 1.8 eV. 
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Figure 1(b) Crystal structure and the charge density evolution of BaFe2Asz2. Plots were generated 
using VESTA |8! (the same below). The isosurface is 0.03 e/bohr*. The yellow color represents an 


increase of the charge density, while the blue a decrease (the same below). œ = 1.8 eV. 
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Figure 1(c) The relationship between the maximum value of charge density change and the 


evolution steps. w = 1.8 eV. 


When w = 0.9 eV for the applied electric field, Figure 2(a) shows the time 
dependence of the electric filed. Figure 2(b) shows the crystal structure and the 
real-time evolution of the change of the charge density. The change of Fe ions 
(marked by 3 and 4) reaches is maximum after 800 steps and the frequency is about 
160 meV. Compared with Figure 1(b), it is found there are significant differences in 


the pattern of charge density changes. 
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Figure 2(a) Time dependence of the electric field applied, œ = 0.9 eV. 
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Figure 2(b) Crystal structure and the charge density evolution of BaFe2As2. The isosurface is 0.02 


e/bohr’. œ = 0.9 eV. 


For the second method to excite the charge vibration, the position of an As atom 
is slightly changed. The original coordinate is (0.0000, 0.0000, 0.3585), and modified 
to (0.0000, 0.0000, 0.3685), ie., the Ca atom was moved in the (001) direction by 0.13 
A. Ground state was calculated. Then the coordinate was modified back to (0.0000, 
0.0000, 0.3585) for TDDFT calculation. Charge densities were recorded after 5000 
and 8000 steps, respectively. Figure 3(a) shows the charge density change recorded 
after 5000 steps’ evolution. The change of the Fe ions marked by an arrow reaches its 


maximum after 500 steps and the frequency is 250meV. Figure 3(b) shows the charge 
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density change recorded after 8000 steps’ evolution. The change of the Fe ions 
marked by an arrow reaches its maximum after 600 steps and the frequency is 
200meV. It can be seen from Figure 3(a) and Figure 3(b) that the pattern of charge 


density changes can be different for the same excitation. 
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Figure 3(a) Charge density change recorded after 5000 steps’ evolution, isosurface =0.05. 
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Figure 3(b) Charge density change recorded after 8000 steps’ evolution, isosurface =0.1. 


The author has also investigated other unconventional superconductors, and 
obtained similar results. For LazFe2As2O2, FeSe sheet and HgBazCa2Cu3Os, the 
characteristic frequencies are 160 meV, 190meV, and 250meV, respectively |1°, 
Orbital (charge density of Fe or Cu ions) fluctuation may also be the electron pairing 


mechanism for copper oxide superconductors. It is worth noting that the change of the 
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electron clouds of transition metal ions is very complicated. There may be different 
modes, corresponding to different frequencies. 

Why do the electron clouds of transition metal ions change this way? The main 
reason is that the 3d shell of transition metal ions is not fully filled. The electron 
clouds possess no spherical symmetry, and easy to change under electric fields. 

The author proposes that the electron pairing mechanism (for both iron-based 
and copper-oxide superconductors) is as follows. When a free electron comes to a 
new place, the electron clouds of the neighboring transition metal ions will change. In 
this way, the charge densities around the free electron will decrease. When the free 
electron leaves, the electron clouds of the transition metal ions will not relax 
immediately (just as the lattice vibration), so that there will be a region lack of 
electron. Another free electron will be attracted. An attraction appears. This 
mechanism is essentially the same as the electron-phonon interaction, except that the 


medium is the change of the electron cloud, not the displacement of the ions. 


4 Summary 

Real-time evolution of the electron clouds of transition metal ions in 
unconventional superconductor was calculated by the TDDFT method. The 
frequencies of evolution of electron clouds match well with the frequencies of the 
lattice vibrations in conventional superconductors at normal and high pressures. 
Though the frequencies obtained by this method are not accurate, it can give a 
significant evidence that the change of electron clouds can be the electron-pairing 


medium in unconventional superconductors. 
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